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Abstract. We present a V light curve of the eclipsing 
binary system HV2543 (Sk-67°117) in the Large Mag- 
ellanic Cloud based on CCD images acquired between 
the years 1995 and 1998. We have analysed this light 
curve and published radial velocity data, finding that this 
system is semidetached, with the secondary (less mas- 
sive and less luminous component) filling its Roche lobe. 
From our analysis with the Wilson-Devinney code, we es- 
timated the following masses and radii for the components 
of HV2543: Ah = 25.63 ±O.7M0, i?i = 15.54 ±O.4R0, 
M2 = 15.63±1.0Mq andi?2 = 13.99±0.4Rq. On the ba- 
sis of B and V photometry of the field stars, we found that 
HV 2543 is member of an OB association, perhaps related 
to which the massive binary system Sk-67°105 belongs. 

Key words: binaries: eclipsing - stars: early-type - stars: 
fundamental parameters - stars: individual: HV 2543 



1. Introduction 

The analysis of light curves of eclipsing binaries, in addi- 
tion to radial velocity data, provides fundamental knowl- 
edge about the masses and physical dimensions of the 
stars. Studying massive binaries in the Magellanic Clouds 
we can learn about the evolution of these systems at metal- 
licities lower than that of our galaxy. 

The Harvard variable HV2543 {a = 5^27^27% S = 
-67°11'54", J2000) is a hot binary star in t he LM C. It 
was catalogued as an OB star by Sanduleak ( 1969 ), who 
assigned it the identification -67° 117. The eclipsi ng na - 
ture of this binary was discovered by Gaposhkin ( 1970|) , 



Table 1. Instrumental configuration 
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Detector 
readout noise 
gain 
filters 

telescope configuration 
focal relation 
scale 
field size 



Tektronix 1024 x 1024 CCD 
10.4 e" 

7.97e-/ADU 

Johnson V and B 

direct CCD -I- focal reducer 

2.83 

0.813 arcsec pixeP^ 
diameter ~ 9', circular 



who published a photographic light curve fitting a pe- 



riod of 4.829052 days (see also Payne-Gaposhkin 1971). 
Photoelectric photometry was performed by Isserstedt 
( |1979D who found V = 12.92, [B - V) = -0.18 and 
[U — B) = —1.03. Radial velocity orbit was obtained by 
Niemela & Bassino ( |1994D who derived physical parame- 
ters of the binary components and concluded that HV 2543 
was a semidetached system with the less massive compo- 
nent filling its equipotential Roche surface. On the ba- 
sis of their spectroscopic data, they classified this system 
as 08V:+09III. Smith Neubig & Bruhweiler (|1999D have 



published an UV spectral classification of LMC OB stars 
based on lUE data, assigning to HV 2543 the type 09III. 

In this paper we present a CCD V light curve for 
HV 2543. By means of the combined analysis of these data 
and previously published radial velocities, we derive new 
values of the fundamental parameters of this system. 

The paper is organized as follows: in Sect. 2 we de- 
scribe the observations, reductions and calibration steps. 
In Sect. 3 we describe the photometric results and light 
curve fitting. In Sect. 4 we discuss the results and in Sect. 
5 we present our conclusions. 
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Table 2. Observing logs 



date 1995, Oct 21-25 1997, Nov 15-21 1998, Dec 2- 

HJD-2450000 0011-0016 0767-0774 1149-1156 

number of images acquired x filter 35 x F 19 x B + 42 x 1/ 36 x V 

range of exp. times [s] 15-40 5-40 7-30 



ences between sky-flats and dome-flats were found to be 
less than 0.3 % at the edges of the frames, and negligible at 
the center. Whenever high quality sky-flats (that is, zero 
clouds) were available, these were preferred to dome-flats. 

For each observing run, the following steps were per- 
formed: 




A master image was made by aligning and combining 
the highest quality frames (best seeing and darkest 
sky). On this combined image, a proflle fitting pho- 
tometry was performed using a stand-alone version of 



DAOPHOT II (Stetson 1987, Stetson 1991). The re- 



sulting output was used as master coordinate list in 
the following steps. 

Aperture photometry was performed on each individ- 
ual frame using a circular diaphragm of 9.5" diameter. 
The sky background was determined locally for each 
star using a sky annulus of radii 9.76 and 20.33". 
A list of instrumental magnitudes for each individual 
frame and for each star in the master list was con- 
structed. 



Fig. 1. Finding chart for HV2543. North is at the top 
and East is at the left. The size of the region displayed is 
approximately 5.4' x 6.4' 



2.1. Observations 

The CCD images here analysed were acquired with the 
2.15-m telescope at CASLEO, during three runs in 1995, 
1997 and 1998. About five frames of the HV2543 held 
were obtained each night, except during eclipses, when 
the observation was more intensive. During each observing 
night, a series of 10 ^ 15 twilight flats and bias frames 
was also obtained. In addition, for the 1997 and 1998 runs, 
10 ~ 15 dome flats were also obtained each night. Tables |] 
and H list details of the observations and Fig. ^ shows a 
flnding chart for IIV2543. 

2.2. Reductions and Photometry 

Reductions were performed using IRAF routines. For each 
night, the frames were reduced in the standard way (over- 
scan and bias corrections, flat flelding) using properly 
combined bias and flat frames. The illumination differ- 



Five faint neighbour stars were detected in the com- 
bined frame within a radius of 9" around IIV2543, being 
the brightest of these some 5 mag fainter than the binary. 
Light contamination from these stars slightly contributes 
to the measured magnitude of HV 2543 (differences of the 
order of 0.014 mag). Consequently, an adequate flux cor- 
rection was applied to the aperture photometry. 

2.3. Photometric calibration 

No single objects were chosen as comparison and check 
stars. Instead, zero point corrections were computed for 
each frame using a group of ~15 stars, covering widely 
the range of magnitudes of HV2543. To put all the obser- 
vations in an unique instrumental system, we proceeded 
as follows: 

1. A synthetic reference system was deflned using the 
mean instrumental magnitudes of the calibration stars. 

2. For each frame, a zero point offset was computed and 
the magnitudes were transformed to the synthetic ref- 
erence system. 

3. After this transformation, the RMS residuals were 
computed for each star. Only the stars with the small- 
est residuals (< 0.01) were conserved as calibration 
stars. 
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4. We repeated these steps with the new, cleaned calibra- 
tion star list. 



to an unified instrumental system. They do not include 
the error of the transformation to the standard system. 



Since all nights were not really photometric, this pro- 
cedure allowed us to have a direct estimate of the internal 
errors of each image. 

The above enumerated steps were performed indepen- 
dently for each observing run. We proceeded in this way 
because the stars available as calibrators arc not the same 
for each run, due to pointing differences. On the other 
hand, some stars remain perfectly stable during a run, but 
show luminosity variations from one year to another, be- 
ing then inadequate as calibrators between different runs. 

After this first stage of calibration was complete, a list 
of robust estimates of instrumental stellar magnitudes and 
their dispersions was constructed for each run. After cross 
identification between the three lists, the most stable stars 
were selected to serve as calibrators between the different 
runs. 

During the 1997 run, 19 frames through the Johnson's 
B band were also acquired. These frames were used to 
determine the transformations to the standard V system. 

The final transformations to the standard system were 
performed by means of aperture photometry of 33 stan- 
dard stars - ranging from (B — V) = —.004 to 2.192 - from 



Landolt ( 1992 ) acquired during the photometric nights. 
We used the following transformation equations: 

b = B + bi+b2 X + b3 {B -V) 

V = V + VI+V2 X + V3 (B-V) 

We adopted the mean values of the extinction coef- 
ficients given by Minniti et al. ( 1989] ) and obtained the 
following values for the other coefficients: 

Nov. 19, 1997: 
bi = 4.2826 ±0.0105 
63 = -0.0208 ±0.0115 
vi = 3.4446 ± 0.0051 
V3 = -0.0456 ±0.0055 

Nov. 20, 1997: 
61 = 4.2986 ±0.0076 
63 = -0.03116 ±0.0086 
vi = 3.4502 ± 0.0057 
V3 = -0.0535 ±0.0067 

The RMS residuals of the V transformations were 
0.017 mag. We first shifted our unified instrumental sys- 
tem to that corresponding to a V frame acquired imme- 
diately before a series of B exposures, on Nov. 19, 1997. 
Then, we applied the standard transformations derived for 
that night. The Nov. 20 transformations were used only 
to check the quality of the first ones. 

Our photometry of HV 2543 is presented in Table ^. In 
successive columns the heliocentric Julian day, standard V 
magnitude, internal errors, seeing and airmass are given. 
The errors (crj) account for the internal photometric errors 
estimated by DAOPHOT and those of the transformation 



3. Results 

3.1. Ephemeris 

We have observed one primary minimum during 1997, at 
HJD= 2 450 768^802 ± 0.001. Combining this value with 
that given by Payne-Gaposchkin (1971) we obtained a new 
ephemerides for HV2543: 

Eo = 2 450 768^802 ±0.001 
P = 4^829 046 ±0.000 004 

We have used this ephemeris in calculating the phases 
for the radial velocity and photometric data. 



3.2. Light and radial velocity curve analysis 

The V light curve of HV 2543 is typical of near-contact bi- 
naries, with different depths of the minima and continuous 
variations of light in out-of-eclipse portions. Also notice- 
able is the O'Connell effect, i.e., the small difference in 
magnitude between both maxima (see Davidge & Milone 
|1984 ). We solved simultaneously the light and radial ve- 
locity curves using the Wilson & Devinney code (hereafter 
WD), which is very well suited for the study of close bi- 
nary systems (Wilson & Devinney 1971, Wilson 1990). 
Unit weight was given to photometric CCD observations 
since they all were collected with the same instrumental 
configuration and have comparable quality. A unit weight 
was also given to most of the radial velocity points except 
for those near eclipses or indicated by Niemela & Bassino 



(1994) as less confident. We also corrected the data points 
corresponding to the phase (p = 0.05, that are inverted 
(i.e., the 08 velocity corresponds to the 09 and vice- versa) 
in Niemela & Bassino (1994). The relation between pho- 
tometric and spectroscopic weights was given through the 
values of the mean standard deviations (sigma) of the data 
which were estimated to be 17 and 12 km s~^ for the ra- 
dial velocities of the primary and secondary components, 
respectively, and 0.02 mag for the light points. Only the 
ratios between these values are significant. 

The radii that result from our preliminary light curve 
analysis indicate that the 08 component is more alike an 



08III star. Chlebowski & Garmany (1991) give a temper- 
ature of 36000 K for an 08III star. Given that Schmidt- 
Kaler (1982) gives a lower temperature for such spectral 
type, we also performed an analysis assigning a 34700 K 
temperature for the primary. Standard values of bolomet- 
ric albedos, A — 1.0 (Rucinski 1969D, and gravity dark- 
ening coefficients, g = 1.0 (Lucy 1976), for radiative en- 
velopes were used. Linear limb-darkening coefficients were 
determined from tables by van Hamme ( |1993| ). These pa- 
rameters were not adjusted. The adjustable parameters in 
our computations were: a (semimajor axis), Vj (systemic 
radial velocity), i (orbital inclination), T2 (temperature 
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Table 3. V light photometry of HV 2543. 
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.156 


0.013 


4, 


.67 


1 
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0.007 


2, 
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1 
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1 
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2, 
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.557 


0.010 
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1 
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15. 
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0.009 
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1 
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0.007 


2, 
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.16 


1 
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Fig. 2. Top: Observed and modelled light curve for HV 2543. Bottom: (0-C) residuals for the light curve. The symbols 
+ , X and O correspond to the data collected during 1995, 1997 and 1998 respectively. 



of the secondary component), q (mass-ratio), fii and 
(modified potential of both components) and Li (lumi- 
nosity of the primary). The first approaches to fit the ob- 
servations were made by means of the light curve (LC) 
program until an acceptable fit to the V light curve was 
obtained. 

Then we proceeded with the differential corrections 
(DC) code. We first adjusted the parameters a, and q 
using only the out-of-eclipses radial velocity data. There- 
upon, we left those parameters fixed and fitted i, T2, 
^2 and Li, using only the light curve. Then we re- 
computed the first set of parameters leaving fixed the 
second set, and iterated this procedure until the correc- 
tions were smaller than their errors. A better solution was 



achieved allowing the fit of q with the photometric data, 
that is, including q within the second parameter set. Mode 
2 (Leung & Wilson 1977 ) which stands for detached sys- 
tems was used at the beginning, but after a few runs of 
the DC program the computations clearly evolved towards 
a semidetached configuration with the secondary compo- 
nent filling its critical Roche lobe. Thus, mode 5 of the WD 
code was employed until the final solutions were found. In 
this mode, the parameter ^2 cannot be adjusted and it 
is set equal to the critical filling-lobe value. The solutions 
never converged to contact configurations. 

Since some systematic differences were obtained be- 
tween the observed and modelled light curves (due to the 
O'Connell effect), we tried to improve the binary mod- 
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Table 4. Model Parameters 



a 


41.5 ± 0.5 Rq 




293.2 ±6kms-i 


i 


89 ± 1° 


q {M2/M1) 


0.61 ±0.05 


Ti 


36000 ~ 34700 K (adopted) 


f2i 


3.363 ± 0.06 


91 


1.00 (adopted) 




1.00 (adopted) 


XI 


0.541 (adopted) 


T2 


29200 ~ 28180 K 


n2 


3.082 ±0.12 




1.00 (adopted) 




1.00 (adopted) 




0.561 (adopted) 


One hot spot solution (in star 1): 


colatitude 


90° 


longitude 


253.4° 


angular radius 


35.8° 


temp, factor 


1.053 


Two hot spot solution: 


For star 1: 




colatitude 


90° 


longitude 


270° 


angular radius 


39.1° 


temp, factor 


1.032 


For star 2: 




colatitude 


90° 


longitude 


0° 


angular radius 


80° 


temp, factor 


1.030 



^ Values for T2 resulting for the above quoted values 

adopted for T\ respectively. 
° x\ and X2 stand for the bolometric limb-darkening 

coefficients. 



elling by including one or two hot spots on the surface 
of the stars. The spots were placed at the equator of the 
stars (colatitude = 90°) while the other spot parameters 
where adjusted with the WD code. The range of the ad- 
justed temperature for the secondary star, T2, depends on 
the adopted value for Ti. To estimate the errors of the 
other parameters we considered the differences between 
the values that arised from different solutions (i.e., using 
the "spectroscopic" value of q, 0.55, or leaving it to in- 
crease until 0.64, to allow a better fit to the photometric 
data). The adopted solution, with q = 0.61, is a compro- 
mise between the optimal fit of the light curve and the 
tolerable deviation from the spectroscopic data. Table ^, 
lists the model parameters. Fig. ^ and ^ show the mod- 
elled light and radial velocity curves, respectively, derived 
from the one spot solution and plotted along with the ob- 
servations, and their corresponding residuals (0-C). 
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Fig. 3. Top: Observed and modelled radial velocity curve 
for HV2543. Bottom: (0-C) residuals for the radial veloc- 
ities. Hollow circles correspond to the primary component 
and filled ones stand for the secondary. 

The final value of the orbital inclination is incidentally 
high, very close to 90°. This fact should produce a total 
eclipse at the secondary minimum, but, as the dimensions 
of both components are quite similar, this is not a notice- 
able feature of the observed light curve. 

Absolute values of the masses, dimensions and bolo- 
metric magnitudes of HV2543 were computed from the 
one spot solution. They are listed in Table ^ where i?i 
and i?2 correspond to mean values of the derived polar, 
back and side radii of the stars. The final configuration of 
HV 2543 shows two components of similar dimensions but 
different masses and temperatures, with the secondary less 
massive star filling its Roche lobe. 

3.3. BV Field photometry 

We combined the B frames acquired during 1997 to make 
a master B image, on which we performed a profile fitting 
photometry. We used this photometry together with that 
performed on the V master frame to generate a colour- 
magnitude diagram, which is shown in Fig. |^. 

In this figure, the presence of an OB association cov- 
ering the whole frame (9' of diameter) is evident. Os- 



trov et al. (1999) found evidence of a stellar association 
surrounding the massive binary system Sk-67°105. The 
relatively small angular distance between the two stellar 
groups (~ 9.5') suggests that they are probably related. 
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Table 5. Star dimensions 



Ml 
Ri 

Mtoi 1 
^og gi[cgs] 

A/2 
R2 

Mbol 2 

log g2[cgs] 



25.63 ±0.7 Mq 
15.54 ± 0.4 Rq 
-9.12 ~ -8.90 " 
3.46 ± 0.04 

15.63 ± 1.0 Mq 
13.99 ± O.4R0 
-7.98 ~ -7.82 ^ 
3.34 ±0.01 



' The range in the derived bolometric magnitudes 
corresponds to the range in the values adopted for Ti. 
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Fig. 4. Color-Magnitude diagram for the field surrounding 
HV2543. The two filled circles represent the magnitude 
estimates for the components of HV 2543. Magnitudes and 
colours are not corrected for extinction. 



We have obtained spectroscopic data for five stars in 
the neighbourhood of Sk-67°105. For these stars, we have 
estimated the reddening by comparison between the in- 
trinsic {B 



V)o colours given by Schmidt-Kaler (1982) 
and FitzGerald ( |1970 ) and the observed colours. We de- 
rived E{B -V) ^ 0.17 ± 0.015, while for HV2543 itself 
we obtained E{B - V) = 0.20. A detailed study of these 
OB associations will be presented in a forthcoming paper. 

4. Discussion 



We note that the radius derived for the 08 component 
is substanti ally l arger than the one obtained by Niemela 
& Bassino ( 1994 ). In fact, the spectroscopic data suggest 
that the luminosity of the 09 star is larger than that of 



the 08 star, which caused Niemela & Bassino to refer to 
the 09 component as "primary" , even though it is the one 
presenting the largest radial velocity amplitude, and con- 
sequently the less massive component of the binary sys- 
tem. Hence, we explored alternative solutions to the light 
curve that could yield a smaller radius for the 08 compo- 
nent. The values i?i - 12.75, R2 ^ 15.31, Mi - 21.91 and 
M2 ~ 19.34, in solar units, provide a good fit to the light 
curve, although they require a value of q near 0.9, which 
is not compatible with the radial velocity data. Hence we 
have discarded this solution, which on the other hand, im- 
plies a rather low value of the distance modulus, namely 
18.18, assuming Ay = 0.62 (see below) and a dopti ng the 
temperature scale of Chlebowski & Garmany ( 1991 ) . This 
value would result even smaller if a lower temperature 
scale is adopted. 

On the other hand, the temperature difference between 
the two components of HV2543 resulting from the light 
curve analysis is larger than that suggested by the corre- 
sponding spectral types. This fact might be due to the dif- 
ference between the spectral features of the unperturbed 
back sides and the heated inner sides of the stars. Also 
should we have in mind that the mass transfer history 
of the system might account for significant departures of 
the He abundances relative to those regarded as normal 
for LMC members, and this effect could influence some 
spectral lines of the secondary component. 

To solve these puzzles, high dispersion spectroscopy 
would be desirable. 

The semimajor axis and star dimensions of this sys- 
tem are alike those determined by Pritchard et al. (1998) 
for HV2241, but the masses are somewhat smaller. We 
presume that these systems have experienced case A mass 
transfer, being now in the slow stage of mass exchange. 
In these cases, the mass gainer should be indistinguish- 
able of a normal star, excepting that it would stand on 
an isochrone corresponding to a shorter age (see Vanbev- 
eren et al. 1998). From the stellar models of Schaerer et 
al. (1993) (for single stars) we found that a 25M0 star 
takes ~ 6.7 Myr until its radius grows to 15Rq, but such 
star would have an effective temperature of only 30000 K. 
The radius and effective temperature derived for the 08 
star are consistent whit those of a single star of some 3.7 
Myrs and ~ 40Mq , completely out of the range of masses 
compatible with the radial velocity data. We note that 
a similar problem arises from the analysis of AB Crucis 
(Lorenz et al. |1994| ). 

It is clear that there are a series of phenomena that we 
do not fully understand, and consequently, our analysis 
can not be considered definitive. The errors given for the 
derived parameters must be considered with caution. True 
errors are not easy to estimate analytically, since they de- 
pend on the importance of phenomena that are not prop- 
erly accounted for, such as wind shocks, uncertainty of the 
adopted temperature scale, radiation pressure effects, etc. 
In fact, the sizes and positions of the spots used to model 
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empirically the O'Connell effect are rather arbitrary, and 
an equally satisfactory solution could be found with other 
parameters, but these details do not affect meaningly the 
derived star dimensions. 

From our photometry we determine (B — V) — — 0.11± 
0.015 for HV2543. This value is somewhat redder than 



—0.18, obtained by Isserstedt (1979). However, given that 
Isserstedt does not detect the variability of HV2543, we 
think that both measures are still in reasonable agreement. 
If we assume a {B - V)q of -0.31 (Schmidt-Kaler |l982|) , 
R = 3.1 (Koornneef |l982|) then it results an Ay = 0.62 for 
HV2543. This value does not depend on which tempera- 
ture scale we adopt, since for the range of temperatures 
of the 0-type stars the {B — V) colours are degenerated. 

Estimating the bolomctric corrections according to 
Massey & Hunter (1998), we derive a distance modulus 
of (m - M)o = 18.31 ± 0.2 to 18.40 ± 0.2, depending on 
the adopted temperature scale. The error accounts for the 
uncertainties in the estimates of R and the bolometric cor- 
rections. This distance modulus must be considered with 
caution, since this system has experienced strong mass 
transfer and exhibits the above mentioned anomalies. 

5. Conclusions 

We have presented a CCD V light curve of the eclipsing 
binary HV2543 (Sk-67°117) in the LMC. The hght curve 
appearance is almost symmetric, with a slight (~ 0.02 V 
magnitudes) O'Connell effect. From BV photometry of 
the surrounding field, we found that HV2543 probably 
belongs to an OB association not previously identified. 

We have analysed this light curve and the published 
radial velocity data using the Wilson-Devinney code, find- 
ing that HV2543 is a semidetached system with the less 
massive and less luminous member filling its Roche lobe. 
The primary minimum occurs when the 08 component is 
behind the 09. From our analysis, we derived fundamental 
parameters for this system, obtaining Mi — 25.63 ± 0.7, 
M2 = 15.63 ± 1.0, Ri = 15.54 ± 0.4 and R2 = 13.99 ± 0.4 
in solar units. Some discrepancies between the present re- 
sults and those derived fro m the previous spectroscopic 
analysis (Niemela & Bassino 1994), should be thoroughly 
addressed in future works. 

We found that this system has experienced significant 
mass exchange, being the present most massive star the 
originally less massive. 
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